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We consider the prompt photon production in pp colhsions using, within the framework of pertur- 
bative QCD, a non-gaussian distribution for the transverse momentum distribution of the partons 
inside the proton. Our description adopts the widely used in the literature factorization of the 
partonic momentum distribution into longitudinal and transverse components. It is argued that 
the non-gaussian distribution of the intrinsic transverse momenta of the partons is dictated by the 
asymptotic freedom as well as the 3D confinement of the partons in the proton. To make this as- 
sociation more transparent we use the MIT bag model, which plainly incorporates both properties 
(asymptotic freedom, confinement), in order to determine in a simplified way the partonic transverse 
momentum distribution. A large set of data from 6 different experiments have been fitted with this 
simple description using as a single free parameter the mean partonic transverse momentum (kr)- 
Surprisingly enough, a perfect fit of the experimental data turns out to require (kr) values which 
are compatible with Heisenberg's uncertainty relation for the proton and decrease almost smoothly 
as a function of the scaled variable z = where pr is the transverse momentum of the final 

photon and ^/s is the beam energy in the center of mass frame. Our analysis indicates that asymp- 
totic freedom and 3D confinement may influence significantly the form of the partonic transverse 
momentum distribution leaving an imprint on the pp ^ y + X cross section. 

PACS numbers: 13.60.Le,13.85.Ni,12.38.Qk 



The production of photons with large transverse mo- 
mentum is an excellent probe of the dynamics in hard 
scattering processes [J, [I]- In particular, the study of 
direct photon production possesses numerous and well 
known advantages, both theoretical and experimental 2- 
0]. In the latter case the main advantage is that pho- 
tons are easier to detect than jets. From the theoreti- 
cal point of view the main advantage is the simplicity 
of the process allowing for an accurate determination of 
the gluon distribution within the proton. In the lowest 
order (0{aas)) only two subprocesses, gq — >■ jq (Comp- 
ton) and qq — >■ jg (annihilation), contribute to high px 
photons. Their characteristic signature is the produc- 
tion of a photon isolated from the hadrons in the event, 
accompanied by a kinematically balancing high-p^ jet 
appearing on the opposite site. In the next-to-leading 
order (NLO) the process associated with the production 
of a photon coming from the coUinear fragmentation of 
a hard parton produced in a short-distance subprocess, 
constitutes a background to the direct photon production 
of the same order in as as the corresponding Born level 
terms [l^ provided that the fragmentation scale is large 
enough. However, the contribution from fragmentation 
remains small (less than 10%) for fixed target experi- 
ments and becomes significant only in inclusive prompt 
photon production at higher collider energies [lOj. Re- 
cently there has been observed a systematic di sag reement 
between theoretical NLO predictions 0, d, HlHlsj and 
experimental data for prompt photon production 

which cannot be globally improved adapting the gluon 



distribution function. Especially for fixed target exper- 
iments NLO approximation shows a significant underes- 
timation of the cross section for some of the measured 
data sets 
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16l | . A similar discrepancy can be observed 
between NLO calculations and the experimental data of 
inclusive single neutral pion production: pp — > tt'^X in 
mostly the same experiments as in the photon case [l8| . 
For the pion production the theoretical description is im- 
proved by taking into account certain large contributions 
to the partonic hard scattering cross section to all orders 
in perturbation theory using the technique of threshold 
resummation ^l^,[l^. The same technique can be applied 
to the photon production by calculating the QCD resum- 
mation contribution to the partonic processes qg — >■ jq 
and qq jg (2^, 21 1. However, the result is a rela- 



tively small enhancement, not enough to compensate for 
the gap between the prompt photon data in fixed tar- 
get experiments ll|, |2l| and theoretical predictions. An 
additional improvement can be achieved by including in 
the theoretical treatment resummation effect to the frag- 
mentation component succeeding in this way a good de- 
scription of UA6 and R806 m data but still failing to 
reproduce the data of E706 [11]. The conclusion of this 
analysis is that resummed theoretical results present a 
residual shortfall in the description of photon and pion 
production data in fixed target experiments. One pos- 
sible explanation of this effect is the existence of a non- 
perturbative contribution associated with intrinsic par- 
tonic transverse momentum kx 0, fill . 12^ . To incorpo- 
rate this effect in the conventional pQCD one assumes 
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a factorization ansatz based on the statistical indepen- 
dence between longitudinal and transverse momenta of 
the partons. In this treatment the distribution of intrin- 
sic transverse momentum gikx) is taken to have a gaus- 



sian form as suggested by the early work of Dalitz [23 1. 
Using the so called fcr-smearing one can fit most of the 
experimental data [l|, 0, • However the method suf- 
fers from two serious disadvantages concerning the val- 
ues of the introduced non-perturbative parameter (kx) 
needed to fit the data: (a) they are incompatible with 
the estimation of the proton transverse radius according 
to Heisenberg's uncertainty relation and (b) they do not 
depend smoothly on physical parameters of the process 
as beam energy or the transverse momentum of the pro- 
duced particle. These shortcomings remain even with 
the inclusion of higher order contributions in perturba- 
tion theory making the /cj^-smearing approach unattrac- 
tive from the theoretical point of view. 



Within the framework of pQCD an alternative scenario 
for the description of pp collisions has been recently pro- 
posed [2^ , incorporating the calculation of the partonic 
transverse momentum distribution using a phenomeno- 
logical quark potential model introduced in the past to 
describe baryonic spectra 27|, |28| . The underlying idea 
in this treatment is to determine the influence of bound 
state effects on the shape of g{kT)- In the calculation 
performed in [26| gikr) turns out to have a non-gaussian 
profile with a characteristic tail. Furthermore, when ap- 
plied to the description of the pp ^ + X process, 
one succeeds to overcome the disadvantages of the gaus- 
sian gi^kx) discussed above. However, there are still un- 
satisfactory issues in this treatment since a connection 
with first principles is missing and the asymptotic free- 
dom property in the partonic dynamics is not taken into 
account. 

In the present work we will argue that a non-gaussian 
transverse momentum distribution with a slowly vanish- 
ing tail, may originate from the fundamental properties 
of asymptotic freedom and 3D confinement of the par- 
tonic degrees of freedom in the proton. In addition, using 
the direct photon production in pp collisions, it will be 
revealed that, within the proposed scenario, the values of 
(kr), needed for the description of all the available exper- 
imental data, are characterized by a remarkable smooth 
dependence on the scaled variable z = ^ where pt is 
the transverse momentum of the produced photon and 
^/s the beam energy in the center of mass frame. The 
results of our analysis clearly support the initial assump- 
tion that part of the residual shortfall in the perturbative 
description of 7 or tt production in pp collisions is asso- 
ciated with the transversal confinement of the partons in 
the protons. 

In the next to leading order (NLO) of perturbation 
theory the differential cross section of the single photon 



production in pp collision can be written as: 



Ej--^{pp J + X) = K{pT, Vs) ^ J dxadxb 



be 



da 



■fa/p{Xa,Q )fb/p{xb,Q ) -j{ab ^ cj)5{s + t + u){l) 

dt 

where f^/p (i = a, 6) are the NLO longitudinal parton 
distribution functions (PDF) for the colliding partons a 
and & as a function of longitudinal momentum fraction 
Xi and factorization scale Q ^ is the cross section 
for the partonic subprocesses as a function of the Mandel- 
stam variables s, f, u [l|. The higher order corrections in 
the partonic subprocesses are effectively included in ([1]) 
through the ii'-factor, appearing in the right hand side, 
which depends on the transverse momentum of the out- 
coming photon and the beam energy [33,]. According to 
our treatment we first attempt to describe experimental 
data using a minimal modification of the NLO pQCD in- 
troducing partonic transverse degrees of freedom through 
the replacement [l], [2§| : 



dxi fiipixij Q 



dxid^kT,ig{kT^i)fi/p{xi,Q'^) (2) 



in the PDF of the colliding partons {i — a,b). To avoid 
singularities in the partonic subprocesses we introduce a 
regularizing parton mass 25|, |29[ with value close to the 
constituent quark mass ruq = 0.3 GeV in the Mandel- 
stam variables occurring in the denominator of the cor- 
responding matrix elements {361] . Contrary to the usual 
treatment where gikx) is taken to be a gaussian, here 
we will use a different form inspired by the MIT bag 
model 3J|. The main advantage of this model is that 



it incorporates the basic features of strongly interacting 
quark matter, namely, asymptotic freedom and 3D con- 
finement in a very simple fashion. Despite this simplic- 
ity, a self-consistent approach to obtain partonic momen- 
tum distributions within the framework of the MIT bag 
model is techniqally very complicated [35| leading to un- 
integrated PDF which are in general incompatible with 
the factorization ansatz ([2]). In addition it is not clear 
how to embed correctly such distribution functions in the 
pQCD scheme ([1]). However, since our approach is purely 
phenomenological, one can use the MIT bag model wave 
function for the ground state of a parton inside the pro- 
ton: 
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to obtain the partonic wave function 0(|fc|) in momentum 



space through the Fourier transformation: 



(fir exp[ifc • r]4'(r) 
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In ([3]) is a normahzation constant, jq, ji are spheri- 
cal Bessel functions, Yooi ^lo, ^ii are spherical harmon- 
ics and E is the energy of the parton's ground state, 
while in ([5]) i? is just the radius parameter of the MIT 
bag model. From 0(|fc|) it is straightforward to obtain a 
rough estimation of the partonic transverse momentum 
distribution through the projection; 

/ + 00 
dk^cbHy/W^fT^mVW^fTkl) 
-oo 

(5) 

The (fc-r) is then promoted to a free parameter by rescal- 
ing of the kr axes, introducing the distribution: 
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with (kx) being the mean transverse momentum corre- 
sponding to gikx) of equation ([5]) and (kx) being the 
fitting parameter and the mean transverse momentum of 
g. The form of the calculated distribution is presented 
in Fig. 1 for a particular choice of {kr)- At this point 
it should be noticed that the slowly vanishing tail of the 
distribution ([S]) , which turns out to be crucial for the suc- 
cessful description of the experimental data, originates 
from the 3D character of the problem, since the involved 
Bessel functions are eigenfunctions of the 3D Laplacian, 
describing the motion of the partons inside the bag. 

Using eas. (|ll2p we describe the data of 6 different ex- 
periments using as a single fitting parameter the (fcr)- 
value. Within the MIT bag model such a variation can be 
justified as a dependence of the bag radius on pr and ^/s. 
In particular we investigate the single photon produc- 
tion data of the 3 CERN experiments R807 (ISR), UA6 
(SPS), NA24 (SPS), as weh as the FNAL experiments 
E704, E706 (Tevatron) and the PHENIX experiment at 
RHIC BNL ^0]. In Fig. 2 we summarize all the con- 
sidered experimental data, shown with symbols, appro- 
priately scaled to improve the presentation. The dashed 
lines represent the theoretical pQCD predictions using 
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FIG. 1: The partonic transverse momentum distribution 
gikr) obtained using the MIT bag model. 



the resummation technique according to [3l| . The best fit 
results using the improved pQCD with the non-gaussian 
intrinsic transverse momentum distribution (Fig. 1) are 
presented by crosses which are not easily distinguishable 
since they coincide with the experimental data. In all the 
performed calculations we use Q = It is clearly seen 
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FIG. 2: The cross section of prompt photon production for six 
different experiments at various energies. The symbols indi- 
cate the experimental data, while the dashed line describes 
the theoretical results using resummation technique. The 
crosses indicate the fit results using the non-gaussian trans- 
verse momentum distribution of Fig. 1. The are not easily 
distinguishable since they coincide with the symbols present- 
ing the experimental data. 

that with a suitable choice of (fcy) an excellent descrip- 
tion of the data is possible. 
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The significant improvement of this non-gaussian 
treatment is that the optimization procedure leads to 
(fcT)-values which form a weh defined, smooth, decreas- 
ing function of the scaled parameter z = ^ . This prop- 
erty is illustrated in Fig. 3 where the best fit values of 
(fcr), for the single gamma production in the experiments 
studied using a non-gaussian /c^-smearing (Fig. 3b), are 
compared with the corresponding results obtained using 
a gaussian distribution (Fig. 3a). It is observed that in 
the non-gaussian case, for 2; — 0, the resulting (fcr)- 
value is compatible with the transverse size of the proton 
according to Heisenberg's uncertainty relation. This is 
in fact expected since in this limit one enters into the 
non-perturbative regime tracing geometrical character- 
istics of the proton through the considered process. In 
addition, for increasing z, a decrease in the contribution 
of the partonic transverse momenta is anticipated, in full 
agreement with the non-gaussian result, since in this re- 
gion the perturbative treatment becomes more and more 
accurate. In contrary, as it can be seen in Fig. 3a, the 
gaussian ansatz leads to absurd behavior incompatible 
with physical intuition. 

At this point it should be noticed that analogous cal- 
culation for ttq production in pp collisions, have been per- 
formed leading to (fcx) values with similar characteristics. 
In a more complete treatment one should also calculate 
the cross sections for DIS, Drell-Yan and pp processes us- 
ing (fcT)-values depending on z as dictated by the analysis 
of the 7 production. 

In conclusion we have proposed a pQCD scheme for the 
description of the existing experimental data of prompt 
gamma production in pp collisions using a non-gaussian 
transverse momentum distribution for the partons inside 
the proton inspired by the MIT bag model. The mean 
transverse momentum of the partons is used as a free 
fitting parameter. It turns out that a perfect fit of all 
existing data is achieved using (A:T)-values which tend 
to become a smooth function of z = where pt is 
the transverse momentum of the outgoing photon and 
•\/s is the corresponding beam energy. In addition the 
range of the required (A:T)-values is in accordance with 
the geometrical characteristics of the proton when non- 
perturbative effects are expected to be visible (z — 0). 
This surprisingly successfuU result suggests that the used 
partonic transverse momentum distribution captures es- 
sential characteristics (asymptotic freedom, 3-D confine- 
ment) of the parton dynamics inside the proton. Thus, 
it looks quite promising the effort to put this purely phe- 
nomenological treatment on a more strict basis within 
the framework of pQCD as well as to extent it in order 
to describe pA and AA collisions by taking into account 
confinement characteristics of the involved partons. How- 
ever these challenging questions are left for a future work. 



1 












□ PHENIX 






□ 


V NA24 




3 




E706 (800 GeV) 
A E706(S30GeV) 




> 

fh 2 
A 


□ 

□ □ 


E704 
+ R807 
* UA6 




V 

1 


A 









□ ++ + + + + +» + V ¥>* * *0« Av* + OA* V * 


a* ■ 



0.00 



0.05 



0.10 0.15 

(a) gaussian 



0.20 



0.25 



?2 



00 0.05 0.10 0.15 0.20 0.25 



+ + + + + +^ +%-0«P7#. «P.o, iw. .OA. V . 



0.00 



0.05 



0.10 



0.15 

V/VT 



0.20 



0.25 



(b) non-gaussian 

FIG. 3: The (kr) used in the description of six different exper- 
imental datasets for single photon production, referred in the 
text, (a) with a gaussian and (b) with a non-gaussian intrin- 
sic transverse momentum distribution. The inset in Fig. 3b 
shows the dependence of (kr) on 2; in a refined scale. 
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